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Abstract 

We calculate the Feynman formula for the harmonic oscillator beyond and at caustics 
by the discrete formulation of path integral. The extension has been made by some authors, 
however, it is not obtained by the method which we consider the most reliable regularization 
of path integral. It is shown that this method leads to the result with, especially at caustics, 
more rigorous derivation than previous. 

1 Introduction 

In optics it is known as the Guoy phase shift [1] that when light waves pass through a focal 
point, the phase jumps -njl discontinuously. This point is called the caustic point, in which 
the intensity of light beam goes to infinity classically. In quantum mechanics this is the point in 
which two or more classical paths join. So phenomena at caustics are observed in many systems 
which are in the similar structure. In path integral formula they make a change in the phase and 
this effect is known as the Maslov correction [2]. 

The harmonic oscillator is one of the most important systems because it is solved exactly, it 
gives the first order approximations of various systems and so on. Path integral formula for the 
harmonic oscillator was obtained by Feynman himself [3]. It is given by 



K(xf,Xi;T)= ^ . -expz-— -\{x] + Xp ) cos coT - IxjXf] 



(1) 



where xj (Xf) is the initial (final) position with the Hamiltonian 

H = + ]-co^x^. 

2m 2 

However it is known that this formula is valid only for a half-period. In every half-period, 
infinite classical paths join at xj or Xf, so the harmonic oscillator is the system in which caustics 
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appear. The extension for any time interval was made by Souriau [4] first. For toT /n {ojT jn], 
it is given by 

K(xf,x,;T) = . —e ^ exp 

\ 2nm\ smo^ri 

and for ojT/n = [toT/n], by 

K(xf, Xj; T) = e-'^'^ 6(xf - (3) 

where [x] means the maximum integer not greater than x. Later Horvathy [5] derived them 
by modifying Feynman's original method and Liang and Morandi [6] by using the eigenfunc- 
tion expansion for the propagator. However there seems to be no derivation by the discrete 
time formulation of path integral (the discrete time method). Therefore we derive the extended 
Feynman formula by this method. 

Path integral is a powerful quantization method and has plain idea, however, in practical 
calculation some difficulties arise, for example, integral measure is not necessarily given. The 
discrete time method is one of the regularizations of path integral and in some cases it is indis- 
pensable [7-9]. The outline of this formulation is as follows. Path integral formula is defined 
by 

KixF,tF;xi,t,) = {xF\e-f-''^''-"^\x,) 
or, if the Hamiltonian has time translation invariance, it is given by 

K(xF,Xj;T) = {xF\e~f'"^\x,) (4) 
where T is the time interval. Then we write as 

K(xf, Xu T) = lim <Xf 1(1 - ^AtHflxi) 

/N-l N 
Y\dxtY\iXj\il-jAtH)\xj^,) 
i=l j=l 

where we have successively inserted the completeness relation 

J dXi\Xi){Xi\ = 1 

and put At = T/N, Xj = x{ jAt), X!\[ — Xp, Xq — X[. In the harmonic oscillator, by making use of 



(^j\Pj) = 




the matrix element is given by 
{xjlil - UtH)\xj.,) 



moj 



2fi sin ojT 



|(x/ + xf^) cos (jjT - IxjXf^ 



(2) 



-I 
-I 



2nTi 



{Xj\pj){pj\Xj-i) 



2nTi 



> exp 



1 ^ 2-2\ 

1 - -Af — + —io Xj 

n ^2m 2 ^ ' 

iA<f£.f.V)^0(A.). 



(5) 
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Then carrying out the p^-integrals, we obtain 



/N-l 
I |Jx,exp[-5j, (6) 

where we have put Axj = xj - xj-i, xj = (xj + Xj-i)/2 and omitted 0(Af) terms to disappear 
in ^ oo. This is the path integral formula for the harmonic oscillator by the discrete time 
method. Naive integration of Q leads to the original Feynman formula In Owe calculate 
© carefully to derive the extended Feynman formula @ and ©. 



2 The Extended Feynman formula 

At first we write © by matrix notation as 



m 



A = 



(2a -J3 

-/3 2a -fi 

-/3 2a 





S = CxAx + 2'bx + c) 

2 At 



^ 




(8) 







, b = -fiXxj, 0, . . . , 0, Xf), c = a{xj^ + x/) 



\ \ 2a -13 

,0 ••• -p 2a, 

, oj^Af ^ , oj^Af 

a=\ — , p=l + —— 

4 4 

where t means transposition. The xrintegrals in ([6]) are the Fresnel integrals and they are eval- 
uated by 



e^"- dx = r 

oo y a —oi 



dx = e a/ — 

a 



{a > 0). 



(9) 



So we need to know numbers of positive eigenvalues and negative one. By easy calculation we 
find the eigenvalues of A are 



or 



Ak = 2ia-/3cos—) (k = 1,2, ■ ■ ■ ,N - I) 



T kn / to At T kn \ 

A = 4 cos^ — + tan-^ — 

" 9A^V A ONI 



2N^ 4 2N> 
without a and yS and the corresponding normalized eigenvectors are 



(.Vk)l 



2 Ik 

— sin— TT {k,l= \,2,...,N - 1) 



(10) 
(11) 

(12) 



where («^.)/ means the / th element of V/,. 

Since cosine decreases monotonously on < x < :7r, we find A/^'s are ordered according to 
k. So if we obtain the zero point Xo(N) of 

fix) = 2{a-/3 cos ^x)eC\0,Nl 

then k less than [xo(N)] corresponds to negative eigenvalues. Since f(x) increases monotonously 
and /(O) < and f(N) > 0, f(x) = has a unique solution 

xo(N) = - cos-' - = — tan-i — (13) 
n j3 n 2N 

where we have used cos"' x = 2 tan"' yjj^- Further by the Maclaurin expansion of arctangent, 
(fT3]) is evaluated by 

^ ^ coT { 1 (oST \i 
,„(A0 = -{!--(-)... 

SO we find 

XoiN) T — . (14) 
n 

Now we can count number of eigenvalues of each sign. Two cases arise whether coT/n is integer 
or not. We put M = [ojT/n] for simplicity. 

(i) ojT/n ^ M 

Because M < ^ < M + 1, for sufficiently large 

M < xo(N) <M+\. 
So M of N - 1 eigenvalues are negative and A^^ - M - 1 are positive. 

(ii) ajT/n = M 

Because Xq{N) < M and Xo(N) 1 M, for large but finite A^, M - 1 of - 1 eigenvalues are 
negative and N - M are positive and in A'^ — > oo, M - 1 remain negative and the Mth goes 
to and A^^ - M - 1 are positive. 

In both cases, for sufficiently large N, At i^O and so det A 0. To deal with both cases together, 
we put 

^^(m (ojT ^ Mn) 
[M - 1 (coT = Mn) ' 

We proceed to calculate ([8]). We put jc^ as the solution of Ax + b = : Ax^ + b = 0. Then 
making use of the translation x = y + Xc, we obtain 

S = —CyAy-M-'b + c). 
2At ' 

Further making use of the orthogonal transformation y = Pz with P = (vi,V2, . . ., v^-i) , we 
obtain 

K{xf,Xi\T) = lim 2exp[^sJ (16) 

N->co \-fl J 
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2= (sis) {n/''-*^**'"'}. ^-^<-*--'*-) 

where we have rewritten z to x. First we calculate Q. Taking account of the difference of the 
Frensel integrals (HI) by signs of eigenvalues, Q is evaluated by 



i 



m 



\2mfiAt 



N/2 



n/ 

k=i ^ 



N-l ^ 

n / 

k=L+i ^ 



e SIS 1-1*1^*- 



m 



InihAt 



(17) 



Making use of the well-known formulas 



r=l r=l 

to ([TT]). we immediately find 

where we have put criN) = (1 Putting CD) into ^7}, we obtain 

Q = 



n-l 



m<x) 



(t\N) - d=^{N) 



2i 



(18) 



Inifi 



V 



o-2(A')-o=^(W) 



2i 



(19) 



Next we calculate S c- We write the determinant of the 'W - 1" dimensional matrix A as 
Div-i. ThenD„'s satisfy 



= 2«D„ -IS'Dn.u Do = 1, Di = 20- 



so we find 



and 



D 



N (t\N) - d=\N) 



N-\ - 



li 



Dn-i 



(Dn-2 fiDN-i 

* * 



* * 

3^-2 pN-3 1 



where we have written ingredients which are unnecessary for later calculations as *. Putting 
these into (fT6l) . we obtain 



mco {o-\N) + d=\N)){x]^ + Xp^) - cr{NW(N)4xjXF 



2n 



cr\N) - d='^(N) 
5 



(20) 



By (fT9l) and (|20l) . the Feynman formula (fT6l) is expressed as 

X exp 



^(Xf, X/; T) = lim 



ma* 



o-2(iV)-o=-(A') 



1 



2^ cr2(A^) - a^iN) 



(21) 



Because cr(N) e''^^'^ {N oo), if coT/2 e N U {0} then cr\N) - 'a{N) (N ^ oo), so we 
need to classify whether cl>T/2 e N U {0} or not. 

(i) ojT/n ^ N U {0} 

In this case no divergence appears. Each factor converges: 

cr^iN) - 'a{N) 2i sin coT, (r^(N) + a^{N) 2 cos (oT, cr{NW(N) ^1 (A^ ^ oo) 
and L = M, so (ED is 



K(xf,xi; T) 



mco 



■e 2 exp 



mo) 



Ifi sino^r 



+ xf^) cos CjOT - IxjXf^ 



(22) 



2mfi\ sino^r 
This is in accordance with 
(ii) ojT/n 6 N U {0} 

In this case L = M - I and now divergence appears. To handle divergence we slightly 
rewrite (l2TI) as 



mo) 



: M-l _ 



K{xF,xf,T) = lim 



X exp 



moj 



—{xr + Xf') - 



m(j) 

f' /- 2 , „ 2 



(xj + Xf - 2z(N)xjXf) 



2n ' l-z\N) 

where we have put ziN) = o^iN) /(t{N) for simplicity. Further we put 



(23) 



u = -^1 ^^(xj + Xf), V = J—-{x! - Xf), 



2n 



2n 



(24) 



then we can rewrite (l24l) to 



K{xf,Xj; T) = lim 



mco e 



/v-oo •>/ nin .Ji^xnW VI 1 +zm y/\l-z{N)\ 



(25) 



For sufficiently large A^^ 



1 Mn Mn l^MnO 1 
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so we can put 

arg (t{N) = ^ - £(.N), £(N) > 0, e(A^) ^ (A^ ^ oo) 

and then we immediately obtain 

z(N) = '^(^^ = (- 1 )^e2,-.(iv) _ ^26) 

By (|26l) we find there are two cases of divergence in (|25]) . 

(a) M e 2N U {0} 
In this case 

1 -z(AO =2sine(A^)e-'(5-^(^» 

and 

z(A^) ^ 1, arg(l - z(A^)) = (A^ ^ oo), 



so the denominator and the numerator of the last factor in (|25l) diverge. 
To handle this divergence we consider not K(xf, xj; T) itself but 



K{xf, xr, T)f{xi) dxi (21) 



where f(xj) is a well-behaved function because path integral formula has essentially 
meaning within integral. The explicit expression of (|27]) is 



lim -,/—-— — exp — , — f(xi)dxj 

N^oo ^| nih ^K^iM L 2 J Vll+z(A^)l Vll-z(A^)l 



c< 

r- [2 e-'^- r ^^(^/(^)) + ^^f^P[-Taffl^^P[-T^l .. .... 

lim -i/^ exp — , f(xj(v))dv 

j-oo V in ^ji^iNW 2 J vii +z(A^)i VI 1 -zm 

(28) 

where we have made a change of the integral variable xi to v. Making use of the 
formula 

CO *J — oo 

with «; e C, I arg wl < n/2, we rewrite (|28]) to 



A'-- J-c V ^^^^ V|l+z(A^)l V|l-z(A^)l 

(30) 

where y in arguments has become ve'^ arg(i-z) gy using one of the expression of the 
5-function 

1 _£ 

lim —=e ' = 6(x), 
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with r = \l - ziN)\, we obtain 

Xco 
e-'^''6{v)f{xi{v))dv (31) 
CO 

where xi{v) is xi{ve~^^) actually, but because of the nature of the 5-function, xi{ve~^^) 
is the same with xi{v). Making a change of variable v to x/, we finally obtain 

Xco 
e'^^^SixF - xi)f(xi) dx, 
CO 

so we find 

K{xf, X,; T) = e-'^''6{xF - Xj) = e-'^^Sixp - (-l)^X/) (32) 

(b) M e N \ 2N 
In this case 

1 +z{N) = 2sine(A^)e-'(5-'^W) 

and 

z{N) ^ -1, arg(l + z{N)) = -| (A^ ^ «')■ 
In the similar way with (HD, we obtain 

K{xf, xi\ T) = e-'^"6{xF + xi) = e'^^^SixF - (-ifxi) (33) 
Putting (|32l) and (1331) together, we obtain 

K(xf,Xi- T) = e-'^"6{xF - {-I)" Xj) (M 6 N U {0}) (34) 
This is in accordance with (O. 

3 Discussion 

In this paper we have derived the extended Feynman formula by the discrete time formulation 
of path integral. As stated in [10], we can observe clearly that when time passes over every 
caustic point, number of negative eigenvalues increases and the phase correction is multiplied 
in the Feynman formula. 

As pointed out in [1 1], if we integrate ^ formally to 

1 _ 1 

VdetA V(-l)^ldetA|' 

we do not know which branch of V-T should be chosen. The Frensel integral seems to be 
indispensable to obtain correct number of eigenvalues in each sign. 

The extensions of the Feynman formula for other systems like a forced harmonic oscillator 
have been made [12-17]. The extensions by the discrete time method will be applicable to these 
systems. 
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